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Abstract 

We explore the extension of chiral perturbation theory to a meson gas out of 
thermal equilibrium. For that purpose, we let the pion decay constant be a time- 
dependent function and work within the Schwinger-Keldysh contour technique. A 
useful connection with curved space-time QFT is established, which allows to con- 

^O ■ sistently renormalise the model. We discuss the applicability of our approach within 

[j a heavy-ion collision environment 

The nonequilibrium dynamics of the chiral phase transition has attracted considerable 
interest during the last few years. One of the original motivations to analyse this regime 
was the suggestion that the so called disoriented chiral condensates (DCC) could form 
^+ during the plasma expansion after a relativistic heavy-ion collision (RHIC), giving rise to 

observable effects such as coherent pion emission ||. Traditionally, this and other similar 
effects have been investigated using O(N) models in which initial thermal equilibrium 
is assumed and nonequilibrium is parametrised as the time dependence of the different 
lagrangian parameters [TJ, |3|, |], |!|. Such time dependence is indeed the only relevant 
nonequilibrium feature when working for instance within Bjorken's initial conditions ap- 
proach, where observables such as the order parameter depend only on proper time in the 
X • central rapidity region ||. 

The use of effective models for QCD is imperative to describe properly the dynamics 
of the plasma in this regime of temperature and energy density. The O(N) models include 
explicitly the a meson and are valid only for two light flavours (Nf = 2). Besides, for 
strongly coupled systems, perturbation theory in these models is not well defined and 
one needs to use resummation methods such as large N. An alternative, not so well 
investigated, is to use the effective chiral lagrangian formalism, whose main advantage 
is to provide a consistent perturbation theory in powers of p/A x - Chiral Perturbation 
Theory (ChPT) [0- where p denotes generically any meson external momenta or field 
derivative, and A x ~ 1 GeV. Besides, it is also valid for Nf = 3. So far, this formalism 
has been applied only in thermal equilibrium, to study the low T (T = (D(p)) meson gas 
and the chiral phase transition || |9|. 

We will review here recent work jnj on the extension of ChPT to a nonequilibrium 



situation. The key idea is to make use of the derivative expansion naturally incorporated 
in ChPT, in order to study the system not far from equilibrium. Of course, in a later 
stage one could apply resummation methods such as large N so as to extend the validity 
of the approach closer to the critical point, as it has already been done in equilibrium ||. 
Our starting point is the nonlinear sigma model (NLSM) where we let the pion decay 
constant-the only relevant parameter to the lowest order in derivatives-be time dependent. 
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As commented above, such time dependence can be thought of as proper time evolution. 
We take the initial time t = 0, which would correspond to a proper time To ~ 1 fm/c, a 
typical hadronisation time in a RHIC environment ||. Thus, we will consider the NLSM 
action 

S[U] = I dt I d 3 x f -^- tr d^(x, t)d»U(x, t) (1) 

Here, C is the Schwinger-Keldysh contour in time, containing an imaginary leg of 
length Pi = 1/Tj, Tj < T c being the temperature for t < 0, where we assume that the 
system is in thermal equilibrium. Thus, /(£ < 0) = f n — 93 Mev to leading order. For 
t > the system departs from equilibrium. Note that, since we choose that departure 
to be instantaneous, f(t) cannot be analytical at t = 0. This feature can give rise to 
discontinuities in the observables, and even extra singularities at t = 0, as it has been 



noted in fTT|. Finally, in the above equation, U(x) is the Goldstone boson field, which 
we can parametrise as customarily in terms of pions for Nf = 2 plus kaons and eta for 



N f = 3 E3- We sna11 restrict here to N f = 2. 

The new ingredient we need to incorporate in the power counting in order to be 
consistent with ChPT is 



/(*) „„(p\ /(*) if(tW 



Pit)- \\ x r p{ty nt) 




(2) 



and so on. Otherwise, we shall keep f(t) arbitrary. One can think of f(t) as an external 
source, to which we will find the nonequilibrium response of the system. Another alter- 
native, which we will not attempt here, is to treat f(t) as a field and solve for f(t) the 
hydrodynamic equations self-consistently. 

Once we have defined our nonequilibrium power counting, we can apply ChPT to 
calculate the time evolution of the different observables. In doing so, we must pay special 
attention to renormalisation. In standard ChPT, one-loop UV divergences coming from 
the 0(p 2 ) lagrangian are canceled by tree level contributions coming from the 0(p 4 ) 
one, and so on for higher order contributions. Such fourth-order action is the most 
general one preserving all the symmetries. On the other hand, it is a well-known feature 
of nonequilibrium field theories that new infinities (time dependent in this case) can 
arise [11], |TI|. Thus, our fourth-order lagrangian must contain necessarily new terms, to 
account for extra divergences. In order to find them, we will make use of a very fruitful 
analogy: the action (TJ) is equivalent to formulate the NLSM on a curved space-time 
background corresponding to a spatially flat Robertson-Walker metric, with scale factor 
a(t) = f{t)/f{0 + ) |TI|. In this way, we can construct the 0(p 4 ) lagrangian in the following 
way: we just raise and lower indices with our RW metric in the standard (equilibirum) 
terms and add those terms coupling U(x) with the scalar curvature R(x) and the Ricci 
tensor R^x) to this order. The latter involve new low-energy constants that are fixed 
by analysing the energy-momentum tensor of QCD at low energies |R| . 

As a first application of our approach, we have analysed the pion decay functions 
(PDF) to next to leading order, i.e, including loops with ([3]) and tree level diagrams from 
the OQo 4 ) lagrangian. In thermal equilibrium, one can define two pion decay constants, 
f%(T) (spatial) and /£(T) (temporal) due to the loss of Lorenz covariance [Q. The same 



happens out of equilibrium, where those constants turn into time-dependent functions. 
After analysing the corresponding Feynman graphs the result can be written as |1(J 



for t > 0, where 



[m? = f(t) [1 + 2f 2 (t) - Mt)} - 2iG (t) 

m} 2 = f 2 (t)[l + f 2 (t)]-2zG (t) 



(3) 
(4) 
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Here, L n and L 12 are the two new low-energy constants that need to be introduced to 
this order [13|. While L\ 2 is already finite, L\\ needs to be renormalised in dimensional 
regularisation. In (|3])-([|), Go(t) is nothing but the equal-time pion two-point function 
G (t) = G (x,x) with G (x,y) the solution of the differential equation 



{n x + m 2 (x )}G Q (x,y) = -6 c (x -y )6®(x-y) 



(6) 



with KMS equilibrium conditions Gq(x, ti — iPf,y) = Gq(x, £«; y), tj < and m 2 {t) = 
—f(t)/f(t) plays the role of a time-dependent mass. It is important to remark that this 
mass is a consequence solely of the nonequilibrium behaviour and has nothing to do with 
a explicit chiral symmetry breaking pion mass term. In fact, our model is exactly chiral 
invariant, since, for simplicity, we have chosen to work in the chiral limit (there are no 
pion mass terms in (|1|)). In the language of curved space-time QFT, m 2 (t) is the minimal 
coupling with the metric preserving chiral invariance. Note that m 2 (t) can be negative, 
thus allowing the existence of unstable long-wavelength modes, which play an essential 
role during the plasma evolution 0, |3|, f|, |fl. 

The results ([5])-([|) reproduce the equilibrium f^T) when we switch off the time deriva- 
tives of fit). We remark that Go(t) contains UV divergences, giving rise to time-dependent 
singularities, to be absorbed by f\{t) and f 2 {t) in the renormalisation of L n . An inter- 
esting consequence of our result is that /*(£) 7^ flit) to one- loop, unlike the equilibrium 
In addition, from ©-(H) and (j|) we see that the difference [f^it]] 2 - [fl(t)\ 2 



m 



case 

is finite, so that we can renormalise both at the same time, which is another consistency 

check. 

Therefore, given /(£), all one has to do to this order is to solve the leading order 
propagator equation (^). Of course, that is not possible in general, as it is well-known in 
the context of curved space-time QFT, where solutions can be analytically found only for 
a very few choices of the scale factor, or equivalently, for f(t) [|l|, [16]]. Hence, in order 
to estimate the nonequilibrium effects within our approach, we shall perform a short- 



time approximation, 
response time and t r 



We will restrict then to + < t < t max where + means a small 
ax ~ l//,r — 2 fm/c. For that range, the expansion in t is also a 



chiral expansion, so that we can truncate it to the relevant order. Thus, solving for the 
propagator by expanding m 2 {t) in power series, we get 



2 



frit) = i 
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H< 



? + 0(p7AJ) (7) 



where T c = y/Ef*, m 2 = -/(0+)//(0+) and H = /(0 + )//(0+). Notice that unstable 
modes appear for m 2 < 0. Here, the two constants /jj are already renormalised in terms 
of Ln and L\i and are fixed by the physical value of /*'*(£ = + ). We remark that the 
t = discontinuities are a consequence of our choice of initial conditions. Nonetheless, 
the effect of the L's turns out to be relatively small since L r u , L\ 2 ~ 10 -3 |T3[ . 



Within our short-time approach, we can estimate the thermalisation time by f n (tf) = 
f (neglecting the effect of the L's) which is also the freezing time, since f w reaches its zero 
temperature value. It is clear that by expanding in t we cannot reproduce a thermalisation 
process where Tf ^ 0. The maximum value tf ~ t max is reached for H < and m 2 < 0. In 
fact, unstable modes always tend to cool down the system. Our estimates are somewhat 
lower that typical O(N) calculations @, ||, which is not surprising, since in those works the 
Ti > T c and then gradients are too high to compare with our approach. For instance, in 
||, Ti ~ 200 MeV and \H\ ~ 400 MeV. Thus, we expect our model to be valid when some 
cooling has already taken place and the system enters the validity range of low-energy 
ChPT. Nevertheless, it is worth remarking that large N methods would allow us to depart 
further from equilibrium. Other extensions and applications of our model to be analysed 
in the future include the long-time evolution - by choosing suitable parametrisations for 
fit)-, the Nf = 3 case and the formation of DCC. The inclusion of quark masses and gauge 
fields would allow to investigate, respectively, the quark condensate time dependence and 
photon production in the pion sector. 
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